Localization-Delocalization Transition in the Random Dimer Model 
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The random-dimer model is probably the most popular model for a one-dimensional disordered 
system where correlations are responsible for delocalization of the wave functions. This is the 
primary model used to justify the insulator-metal transition in conducting polymers and in DNA. 
However, for such systems, the localization-delocalization regimes have only been observed by deeply 
modifying the system itself, including the correlation function of the disordered potential. In this 
article, we propose to use an ultracold atomic mixture to cross the transition simply by externally 
tuning the interspecies interactions, and without modifying the impurity correlations. 

PACS numbers: 64.60.Cn,03.75.-b,67.60.Bc 



In a one-dimensional disordered system, Anderson lo- 
calization is known to occur at any energy when the 
disorder is i5 correlated [H, Nevertheless, if one in- 
troduces particular short-range correlations, delocaliza- 
tion of a significant subset of the eigenstates can appear. 
This happens in the random-dimer model (RDM) [3], in 
which the sites of a lattice are assigned energies Ca or ef, 
at random, with the additional constraint that sites of 
energy ef, always appear in pairs, or dimers. The same 
occurs in its dual counterpart (DRDM) 3] , in which lat- 
tice sites with energy ei, never appear as neighbors. In 
these models, extended states arise from resonant modes 
of the (dual) dimers which present vanishing backscatter- 
ing at energy i^rcs- In the thermodynamic limit, the ratio 
^/N /N between the number of delocalized states and the 
total number of states vanishes, and there is no mobil- 
ity edge separating extended and localized energy eigen- 
states. Nevertheless in finite size systems, thus in real 
systems, a localization-delocalization transition can be 
induced by driving i^ios inside the spectrum. 

This model was proposed to be the possible mecha- 
nism which leads to the insulator-metal transition in a 
wide class of conducting polymers such as polyaniline 
and heavily doped polyacetylene (see, for instance, (Jj) 
and in some biopolymers such as DNA Q . The evi- 
dence of delocalized electronic states was experimentally 
demonstrated in a random-dimer GaAs-AlGaAs super- 
lattice 01 , while for photons, a RDM dielectric system 
was used [8[. Recently, a RDM setup has been proposed 
to demonstrate the delocalization of acoustic waves 
For polymers, semiconductor lattices, photonic crystals 
and elastic chains, the dimer resonant energies cannot 
be modified without changing the sample itself. Thus 
the localization-delocalization transition for a (D)RDM 
chain as a function of the relative position of the resonant 
modes with respect to the band modes cannot easily be 
studied using these physical systems. 

In this article, we propose an experimental procedure 
to realize a DRDM experiment with a one-dimensional 
(ID) two-component ultracold atomic mixture in an op- 



tical lattice, and we demonstrate that the localization- 
delocalization transition can be explored by tuning the 
interparticle interactions. 

To introduce disorder, a component (Bd) has to be 
classically trapped in the minima of the potential 
12l |. For this purpose, one can choose a spin-polarized 
Fermi component or a strongly repulsive hardcore Bose 
gas. The other component (J5/) must be able to tunnel 
through the potential maxima. A single impurity Bd 
trapped in a lattice site causes an energy shift of the 
effective potential seen by the second species Bf with 
respect to the case where the impurity is absent. In the 
following we will focus on a boson-boson mixture, taking 
recent experiments on the ^^K-*''Rb mixture at LENS 
as a guide [l3l |. This mixture has tunable interspecies 
interactions for both ^^Rb and in the \F = l,mf = 



1) state 1J|. The ^^K condensed component plays the 
role of the "tunnelling bosons" Bf and the heavier ^^Rb 
atoms are the defects Bd- 

To study the effect of correlated impurities Bd on 
matter-wave transport, we use the ID effective tight- 
binding (TB) Hamiltonian for bosons Bf, 

Us ris — l 

where Ug is the number of sites, Ei G {ea,eb}, and eb's 
never appear as neighbors. The hopping term ti can take 
the values taa between two sites with energy ea or tab be- 
tween two sites with different energies. The constraint on 
sites with energy eb fixes the hopping energies tab to be 
distributed as "dual dimers" of the form Ca £& Ca- 
To realize such an experiment, one needs to assure that 
there can be either zero or one impurity Bd in each lat- 
tice site and that they never appear in succession. The 
procedure we propose is schematized in Fig. [1] (i) First, 
Nb^ atoms of species Bd are trapped in a lattice with 
a step 3d, with the condition that no site must be dou- 
bly occupied, as in the Tonks gas experiment described 
in jTsj. Such a lattice can be realized using beams of 
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FIG. 1; (Color online) Experimental scheme to generate the 
correlated disorder, (a) Impurities are trapped in a lattice of 
step 3d. (b) A second lattice of step d is switched on, and 
the first one turned off. In the final configuration (c), the 
impurities follow a DRDM distribution. 



wavelength A = 2d, tilted by an angle a = acos(l/3) 
[see Fig. [Ija)]. For the case A = 800 nm and a lattice 
potential depth Ug"' ~ 30Eb^, Eb^ = An'^h^ / {2mB^\^) 
being the recoil energy for a boson of mass rn_Bd , and 
for an axial confinement of 2-k x 60 Hz [ ISJ . double oc- 
cupancy can be avoided if Nb^ < 20 [16[. The number 
of impurities can be increased by increasing the lattice 
potential depth or by relaxing the axial confinement, (ii) 
Then impurities are forbidden to occupy neighbor sites 
by adiabatically ramping up the power in a second su- 
perimposed lattice of step d, and switching off the first 



one [17[. The potential depth Ug^ of the final lattice 
UBd,Bfiz) — sin^(7rz/d) must be large compared 

to Eb^ to neglect the impurity mobility during the ex- 
periment course (0.5 to 1 s [3, This condition can 
be fulfilled at Ub^ — ^SEb^ in the presence of attractive 
interactions with the species Bf. Differently, the depth 
U'^^ for the species Bf must be > 2 EBf, the recoil en- 
ergy for a, Bf boson, to guarantee the validity of the TB 
description. 

The effective Hamiltonian (jlj is obtained by a ID 
reduction of the system Hamiltonian by introducing 
the transverse widths (yA_Bf,Ba of condensate and 
of the impurities wave functions in a cylindrical trap 
(20j . Using a TB scheme we introduce the Wan- 
nier function 4'i{z) approximated by the Gaussian func- 

tion0,(z) = [<|^m/{^y^al%)]eM~{^~^^?|{'2^lB,)]. 
where |0i(O)p is the number of bosons Bf \n the lat- 
tice well i. Similarly, the density of impurities is ns^ oc 
exp[— (z — Zi'Y/alB ]■ The determination of the 
widths (J±Bj,Bd a-nd (JzBj,Bd is carried out variation- 
ally 0,0-' 

We can now evaluate the parameters entering the ef- 
fective Hamiltonian ([T|). The site energies are given by 



E,, 



dz 4>iiz) 



2mBf 



+ Ub;{z) 



(2) 



(i>t{z) 



h'^/(2mBfO-'^Bf) + h^Bfi^]_Bf^lBf^ andw^Bj. is the ra- 
dial frequency of the harmonic trapping potential. 4'i{z) 
are modified Gaussian functions, obtained by imposing 
the condition J (f)i{z)(f)j{z) — Sij [21]. The parameters g 
and g' are the strengths of the ID BfBf and BfBj^ inter- 
actions, which are given hy g — (Airh'^a) / {2iTmBfCr'^ Bf) 
and g' — {2'Kh^a') /['Kmr{(j\ b^- with a, a' the 

BfBf and the BfBd scattering lengths and rur the BfBd 
reduced mass. The hopping energies ti are given by 



ti = / dz(t)i{z) 



2m f 



+ UBf{z) 



(i)i+i{z). (3) 



This completes the determination of the effective ID 
Hamiltonian for bosons Bf. 

Delocalization occurs for energy values near to the 
resonance energy of a single dual dimer embedded in 
a perfect lattice of site energies and hopping ener- 
gies taa- For such Hamiltonian H, the wave function 
at energy i? is \ip) = |fc) -f G°T|fc), where \k) is the 
wave function of the unperturbed periodic Hamiltonian 
^^0 = Y.n=-o.^a\n){n\ + taa{\n){n + C.C.), G° the 
unperturbed Green's function G'^{E) = {E — Hq)~^ , and 
T the matrix T{E) = Hi{l - G°Hi)-^, where Hi is a re- 
mainder defined a.s Hj = H — Hq — {eb~ea)\0){0\ + {tab — 
taa){\ — 1)(0| + |0)(1| +C.C.). Here and below the complex 
energy E is considered in the limit of vanishing posi- 
tive imaginary part. Using the renormalization scheme 
outlined in [22] the scattering T matrix in the subspace 
{] - 1), 11)} can be written as T = Hiil - G°Hi)-^, Hi 
being the renormalized remainder Hamiltonian. We find 
Hi = a{\\) with a = tlJ{E - e^) - tlJ{E - e,). Thus, 
the scattering matrix on the subspace {] — 1), ]!)} is iden- 
tically null if a = 0. This occurs at the resonance energy 



, j-2 _ , +2 



(4) 



where ruBf is the mass of the boson Bf, GBf 



Eigenstates are delocalized if i^rcs is inside the lowest- 
energy band E{k) — ea + 2taaC0s{kd), namely if [Ae] < 
2\t'^ - 1], with Ae = (ef, - ea)/taa and t = tab/taa, as 
found by Dunlap and collaborators The correspond- 
ing phase diagram in the {t, Ae) plane is shown in Fig. [2] 
The central lobe corresponds to the case where the pres- 
ence of an impurity is disadvantageous to the hopping of 
a boson Bf. In this region delocalization occurs only for 
values of eb near the center of the band. The left side 
of the diagram corresponds to an increase of the hop- 
ping probability due to the impurities. In this region 
the energy value tb can be in the gap; thus, the disorder 
strength can be very large, but delocalization is estab- 
lished just by the value of i?ros- If -E'rcs is not an energy 
value of the spectrum, for a sufficiently long lattice, all 
states are localized. 

For the evaluation of the site and hopping energies, 
we consider a system of Nb^ = 1.3 x 10'' ^"'K atoms 
distributed in 200 wells, 10% of which are occupied by 
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FIG. 2: (Color online) Phase diagram of the DRDM in the 
plane (t, Ae), with t = tab/taa and Ae = (et — ea)/taa and 
calculated trajectories for different number of atoms Nsf 
obtained by varying the interspecies scattering length a'. 
Dashed line, Nbj = 1-6 x 10*, a £ [0, — 87ao]; continu- 
ous line, Nbj = 1.3 x lO", a £ [0, -380ao]; dotted line, 
Nb; = 1.0 X 10'', a e [0, -570ao]. With 1.3 x 10* atoms, 
the localization-delocalization transition for the BfBd mix- 
ture occurs for a = — 346ao (point C). Other points corre- 
spond to a = -316 (A), -331 (B), -361 (D), and -376 ao 
(E), respectively, ao being the Bohr radius. 



a ®^Rb atom. We choose the depth Ug^ equal to 2.5 
and the optical lattice wavelength A — 800 nm 
(red detuned for both species). For linearly polarized 
beams, this fixes the potential depth Ug for the defects 
to 18 Eb^, and impurity tunneling time of the order 
of 1 s for points A to E in Fig. The effect of the 
interactions is enhanced by a tight radial confinement 
i^_LB//27r = 60 kHz. For such a system the phase dia- 
gram can be explored just by varying the BfBd scattering 
length a' (lines in Fig. [2]). In the experiments this can 
be done by exploiting interspecies Feshbach resonances 
14| . The point (1,0) corresponds to a' = 0: species Bf 
do not interact with impurities; thus, neither site nor 
hopping energies are modified by the presence of Bd and 
the lattice is not disordered. Higher points of the curve 
correspond to greater and greater attractive BfBd inter- 
actions. 

We study the spectrum properties across the transi- 
tion by evaluating the Lyapunov coefficient ^{E), which 
is equal to the inverse of the localization length £{E), 
through the asymptotic relation 



7{E) = m)] 

where G{E) 



lim In 

is-!-oc nsd 



Gn 



XE) 



Gl,nAE) 



(5) 



{E - He 



is the Green's func- 



tion related to the Hamiltonian Hbj at energy E, and 
G^^J{E) = {i\G{E)\j). The matrix element Gi^nAE) 
and Gn^^Us have been computed by exploiting a renor- 
malization/decimation scheme |23j. The behaviour of 



in Fig. [3] The different lines, which correspond to the 
crosses in Fig. [2l show that the localization length is 
greater than the system size for points C, D, E. The lo- 
cation of the minima corresponds to the position of the 
resonance energy -Ercs, which moves inside the band for 
increasing values of \a'\. The nonzero value of 7 in the 
delocalization regime is due to the finite value of in 
computing Eq. ^ (ris=1000 for the evaluation of 7). 

The nature of the states determines the matter-wave 
transport properties near equilibrium. These can be eval- 
uated by embedding the whole system in a perfect lattice 
of site energies ea and hopping energies taa [22|, as pre- 



viously outlined for the evaluation of the dual-dimer res- 
onance energy. The transmission probability T = |tP, 
defined as the squared modulus of the transmission am- 
plitude, 



^ — 2zfc(ns — l)a 



T 



+ iTi^i, 



(6) 



the Lyapunov coefficient through the transition is shown 



for different values of a' is shown in the inset of Fig. [31 
For a' — —346 qq the resonance fits in the band-edge, and 
the corresponding transmission peak arises. For a' = 
— 361ao and a' — — 376ao the peak moves toward the 
center of the band in agreement with the position of the 
minimum value of 7. The width of the peak decreases 
by increasing the system size, as the percentage of the 
delocahzed states scales as i/nj/us. 

Since the condensate energy corresponds to the low- 
est allowed energy (quasimomentum fc = 0), the region 
nearby the resonance could be explored by preparing 
Bloch states with initial quasimomentum 7^ by intro- 
ducing a constant frequency shift between the two waves 
generating the lattice [2J]. We expect that, in the lo- 
calization regime, for any k G [—Tr/d,Tr/d], the whole 
condensate stays at rest in the reference frame of the 
moving lattice, while, in the delocalization regime, at 
k = k^Ej-ce), the bulk of the condensate stays at rest in 
the laboratory reference frame. 

Both the Lyapunov coefficient and the transmittivity 
show some small peaks (sinkings). These structures are 
due to an underlying order present in the procedure il- 
lustrated in Fig. [T] In fact, even if our proposition al- 
lows the distance between two subsequent impurities to 
be equal to any integer > 1, still every three sites is 
definitly without an impurity. The evidence that this 
underlying order does not affect the DRDM physical ef- 
fect is shown in Fig. 21 The transmittivity peak for the 
DRDM pattern proposed in this article (Fig. [Ij is in cor- 
respondence with the transmittivity peak for a genuine 
DRDM. Such a disorder pattern could be generated by 
using a dipolar gas for which repulsive interactions may 
avoid next- neighboring occupation 25|. However, at the 
moment of writing, no dipolar gases have yet been cooled 
down to the degenerate regime in mixtures. 

For completeness of our analysis, we compare the two 
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FIG. 3: (Color online) Lyapunov coefficient 7 in units of l/L, 
L — Usa being the lattice length, as a function of the energy 
of bosons Bf. The dot-dashed green line corresponds to point 
A in Fig. [2I the short-dashed blue line to B, the long-dashed 
magenta line to C, the dotted black line to D, and the continu- 
ous red line to E. The inset shows the corresponding behavior 
of the transmittivity T. 
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FIG. 4: (Color online) Transmission coefficient T for TV = 
1.3 X 10* and a' = — 361ao, and different disorder patterns: 
a genuine DRDM lattice (dashed blue line), a DRDM lattice 
generated with the procedure illustrated in Fig. [T] (continuous 
red line), and an uncorrelated lattice (dotted black line). The 
inset shows the corresponding density of states. 



DRDM models with a lattice where the position of impu- 
rities eb are uncorrelated (26j . In this case the transmit- 
tivity drops (Fig. \^ and becomes vanishing for longer 
chains. The residual peak is a signature of the pres- 
ence of a few dual dimers, and it would be washed out 
in an ordinary disorder model with uncorrelated on-site 
and hopping energies. It is worth pointing out that the 
impurity distribution deeply modifies the nature of the 
states, but not the spectrum itself at this low impu- 
rity concentration. Indeed, the density of states (DOS) 
is essentially the same for the three cases as shown in 
the inset of Fig. 2] When the percentage of impurities 
is increased, the underlying periodicity, which is differ- 
ent in the three models, leads to fragmentation of the 
DOS in three, two, or one band. The DOS, U{E), has 
been evaluated by using the Kirkman-Pendry relation 
N{E) = ^lm{[d\nG,^nAE)]l{dE)} 0. 

In conclusion, in this work we show that, at fixed corre- 
lation function among defects, localization-delocalization 
can be induced by varying the impurity cross section. In 
an ultracold boson-boson mixture, where one component 
plays the role of correlated impurities, the rule being that 
no next- neighboring impurities are allowed (DRDM), this 
can be realized by driving the interspecies interaction by 
means of Feshbach resonances. This is a unique oppor- 
tunity compared to other physical domains where this 
class of disorder was previously identified as the possi- 
ble explanation of the mechanism causing the amazing 
conducting properties of disordered ID systems, such as 
conjugated polymers or biopolymers. 
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